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We show that the ideal relativistic spinning gas at complete thermodynamical equilibrium is a 
fluid with a non- vanishing spin density tensor a^. After having obtained the expression of the local 
spin-dependent phase space density /(x, p) aT in the Boltzmann approximation, we derive the spin 
density tensor and show that it is proportional to the acceleration tensor fi M „ constructed with the 
Frenet-Serret tetrad. We recover the proper generalization of the fundamental thermodynamical 
relation, involving an additional term — (l/2)f2 M „<j M '' / . We also show that the spin density tensor 
has a non-vanishing projection onto the four-velocity field, i.e. t^ — a^utt" 7^ 0, in contrast to 
the common assumption t M = 0, known as Frenkel condition, in the thus-far proposed theories of 
relativistic fluids with spin. We briefly address the viewpoint of the accelerated observer and inertial 
spin effects. 

I. INTRODUCTION 

The study of relativistic fluids with spin has a long history. Among the early major contributors to the subject are 
Mathisson [1] and Weyssenhoff [2]. The latter formulated a model enforcing the condition of vanishing projection of 
the spin density tensor over the four-velocity field, namely: 

= o-^u" = (1) 

also known as Frenkel condition. This auxiliary constraint onto the fluid motion did not seem justified by established 
physical principles and for a system of pointlike relativistic massive particles in general does not apply [3]. This 
condition was in fact critically reexamined in a later work by Bohm and Vigier [4] where they proposed to extend 
the theory to a more general motion allowing the four-vector t ^ 0, yet without obtaining a conclusive quantitative 
formulation. A major step in formulating a theory of relativistic fluid with spin was made in 1960 by F. Halbwachs [5] 
who worked out a Lagrangian variational approach, yet keeping the Frenkel condition as auxiliary equation of motion. 

The hydrodynamics of relativistic fluids with spin has since become the subject of specialized literature in General 
Relativity [6], especially in the framework of Einstein-Cartan gravity theory [7], Most of these studies assume the 
Frenkel condition (1), which implies a considerable simplification in the calculations. Since this condition, as has been 
mentioned, does not derive from any major physical principle, it would be interesting to test it in a simple situation. 

Indeed, a system which allows to perform a theoretical independent test exists: the macroscopic ideal spinning 
gas of particles with spin at full thermodynamical equilibrium. The condition of full thermodynamical equilibrium 
implies that this study can be carried out within the framework of statistical mechanics, independently of any previous 
formulation of a hydrodynamical theory and freely from the danger of begging the question. 

One of us studied in a previous work [8] the physics of the ideal rotating relativistic gas of particles with spin (in 
the Boltzmann statistics) at full thermodynamical equilibrium obtaining the expression of the polarization four- vector 
and momentum spectra. In this paper, we resume that study and derive the local expression of all thermodynamical 
quantities for this system seen as a fluid. Particularly, we will obtain the expression of the spin tensor and the spin 
density tensor and show that the latter enters in the general thermodynamical equation relating proper entropy density, 
proper energy density and pressure, coupled to the acceleration tensor. The nature of this additional term, which has 
been already extensively used in the aforementioned works [6, 7], is that of a quantum correction to thermodynamics 
owing to the presence of spin as a further degree of freedom of particles. We will provide the expression of the 
single-particle spin-dependent phase space density showing that it is a non-diagonal matrix in spin coordinates. We 
will finally demonstrate that even for this simple system (indeed the simplest among all possible relativistic fluids) 
the four-vector t defined in (1) is non-vanishing, thus violating the Frenkel condition. 

The paper is organized as follows: in Sect. II we summarize and further explain the main results obtained in the 
work [8]. In Sect. Ill we calculate the local thermodynamical densities of the rotating gas viewed as a macroscopic 
fluid. In Sect. IV we introduce and derive the angular velocity tensor for a rotating fluid and in Sect. V the spin tensor 
showing its proportionality to the acceleration tensor. In Sect. VI we briefly take the viewpoint of the accelerated 
fluid cell in the rotating gas. 

Throughout this work, we will confine our attention to a special relativistic fluid, i.e. we will disregard gravitational 
effects and assume a flat Minkowski spacetime. 
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Notation 



In this paper we adopt the natural units, with K = c = K = 1. 
The Minkowskian metric tensor is diag(l, — 1, — 1, —1); for the Levi-Civita symbol we use the convention e 0123 = 1. 
We will use the relativistic notation with repeated greek indices assumed to be saturated. Contractions will also be 
denoted by a dot for both vectors and tensors; the saturated index in a contraction between a vector and a rank 2 
tensor is specified by the position of the vector, either on the left or right. Double contraction of rank 2 tensors A 
and B will also be denoted by a colon, i.e. A : B. 

Space-time linear transformations (translations, rotations, boosts) and SL(2,C) transformations are written in serif 
font, e.g. R, L; [p] will denote the Lorentz or SL(2,C) transformation taking the timelike vector rat = (m, 0,0,0) to 
four-momentum p. Operators in Hilbert space will be denoted by an upper hat, e.g. R. Unit vectors will be denoted 
with a smaller hat e.g. p; even though the notation is unambiguous, we will make explicit mention of either possibility 
whenever confusion may arise. 



II. ROTATING RELATIVISTIC GAS AT EQUILIBRIUM 

A macroscopic system with given angular momentum at complete thermodynamical equilibrium can only have a 
rigid velocity field: 

v = u x x (2) 

where u> is a constant vector, as seen from the inertial system where the system is at rest. This is a simple consequence 
of enforcing statistical equilibrium and is a well known result for non-relativistic systems [9] which can be easily 
extended to relativistic ones with the additional bound x x| < 1 (see Appendix B in ref. [8]). Obviously, 
for an isolated system, u> can be constant only if it is directed along one of the system's principal axes. Hence, 
thermodynamical equilibrium of a rotating body requires u) to be parallel to J and, therefore, that the system's shape 
is symmetric for rotation around oj. 

The whole physics of an equilibrated system with finite angular momentum J, chemical potential p, and temperature 
T is fully contained in the density operator pj: 

p J = -L cxp [(-H + pQ)/T}PjP v (3) 
Zj 

where Pj is the quantum Wigncr projector onto a given state \J,M): 

Pj = W(rsf dgD J (g-i)%g (4) 
Jsu(2) 

Zj is the partition function: 

Zj =trexp[(-# + pQ)/T]P 3 P v (5) 
and Py is the projector onto localized states \hy) 

P v = Y,\hv){hv\ (6) 

which form a complet set of quantum states for the system in its finite region V. In a first quantization scheme, these 
states are a complete set of functions which vanish outside the domain V. On the other hand, in a quantum field 
theory, there are some complications which have been discussed extensively in refs. [10, 11], yet they do not essentially 
modify the first-quantization scenario (see also Appendix A). 

The density operator (3) defines a grand-canonical ensemble with fixed angular momentum. If both V and J are 
large, it can be shown [8, 11] that this ensemble is equivalent (at least for first moments of statistical distributions) 
to the rotational grand- canonical ensemble whose density operator reads: 

p 0J = ^exp[(-H + p J Q + u-J)/T]P v (7) 

being u> the angular velocity vector and Z^ the relevant partition function: 

Z^ = trexp[(-# + pQ + u ■ 3)/T]P v (8) 
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The transition from (3) to (7), or from (5) to (8) can be carried out through a large J limit of the Wigner rotation 
matrices in (4) and a subsequent saddle-point expansion in the angular variables defining the SU(2) parameters in 
the integral (4); for details see refs. [8, 11]. For studying macroscopic rotating systems, the ensemble (7) is thus the 
most natural choice. It can be shown that indeed, as long as co/T <C 1, the partition function (8) can be written as 
the product of grand-canonical partition functions of infinitesimal cells moving with a velocity (2), according to the 
general wisdom that equilibrated macroscopic systems must be rigidly rotating. 

It is worth showing this in some detail for an ideal relativistic gas of massive particles with spin S in the limit 
of Boltzmann statistics, i.e. considering particle as distinguishable objects. In this case, we can simply write the 
partition function (8) as a product of single-particle partition functions and the calculation of (8) in fact reduces to 
compute matrix elements of single-particle compatible operators h,q,j, like this: 

(p,r\ exp[-(h + fiq + w ■ ])/T]P v \p, a) = C - e ' T+ ^ T {p, t\ exp[u; ■ J/T]P v \p, a) (9) 

where a and r label polarization states . The above equality requires that, as mentioned, [ft,, j] = [q, j] = [q, h] = 0. 

The matrix element on the right hand side of eq. (9) can be calculated by means of an analytical prolongation 
from imaginary co values. After replacing lj/T with — i<f> we have a rotation R&((f>) around the axis u> of an angle <ft 
appearing in the eq. (9). We can then expand the matrix element on the right hand side as: 

(P,t\R&(4>)Pv\p,v) =J2 [ d V (P,r\%m P \a')(p,<j'\P v \p,<j) (10) 

Note that this expansion would have been impossible for the operator exp[u> • j/T] itself because, for real u>, all its 
matrix elements calculated with momentum eigenstates turn out to be infinite. In other words, the correct procedure 
to make the analytic continuation of the matrix element is, first, go to imaginary u> values and thereafter insert an 
identity resolution like in (10). The matrix element of the representation of a rotation is well known; it involves a 
Dirac's delta and a Wigner rotation matrix: 

(p,r|R^(0)b^^-5 3 (p-R^(0)(p')) J D s ([R^(^)(p'r 1 R^(0)[p']w (ii) 

where [p] is defined as a Lorentz transformation (in fact its SL(2,C) correspondant) taking the timclike vector mt = 
(m, 0, 0, 0) to p and D s denotes, as usual, the (2S+l)-dimensional SL(2,C) representation (what is commonly labelled 
as the (S,0) in fact). Obtaining the matrix element of the projector Py over momentum eigenstates is not so 
straightforward. The calculation can be done in a quantum field theoretical framework [11]: 

(p',a'\P v \p,a) = \ ^1 ^d 3 xc-(P-P'> + ^(b'lW^U (0|Pv|0> (12) 

where the last factor (0|Py|0), the vacuum expectation value of the projector Py, is an immaterial factor [10, 11] and 
tends to 1 for large volumes V as Py — > I; we will henceforth neglect it. The formula (12) gives the correct infinite 
volume limit Py — > I, that is the momentum eigenstate normalization. For finite volumes, a striking feature is the 
presence of non-trivial matrix elements in spin coordinates; their particular form can be justified also by requiring Py 
to commute with rotation operators (see Appendix A). 

Plugging (11) and (12) into (10), and neglecting Py vacuum expectation value, one obtains: 

(p,r\R^)P v \p,<j) = jfd'xe*^'*'" 1 ^ ^ (^ S (b]"%Wb]) + £>S (b] tR <i(^)b] t_1 )) Tff (13) 
In order to get to this expression, advantage needs to be taken of the unitarity of the Wigner rotation, i.e.: 

^([R^XpOr^^b'D - ^([MW^M^bf- 1 ) (14) 

and the unitarity of R itself as an SL(2,C) matrix. The equation (13) can now be analytically prolonged to imaginary 
angles yielding the final expression for the matrix element in (9): 

(p,T|exp[ W -j/T]Pyb,<7> = J d 3 xc-(P- R -(-/ T )"(P» i(^(b]- 1 R^(z W /T)b]) + J D s ([p]tR 1 ,( JW /T)b]t- 1 )) Ta 

(15) 



For the spin quantum numbers a we follow the convention in [12], where they are defined as eigenvalues of the projection of the Pauli- 
Lubanski vector divided by the mass onto the third space-like four-vector orthogonal to the particle four-momentum. For more details 
see ref. [12] and [11]. 
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The matrix element (15) is the key ingredient to calculate the equilibrium single-particle phase space distribution and 
thence, the thermodynamics of the rotating Boltzmann gas; it is then worth some comment. 

Firstly, it should be noted that if the system's shape is symmetric by rotation around the u> axis, as has been 
mentioned, the projector Py commutes with the corresponding generator in the Hilbert space: 

[Pv,j-w]=0 (16) 

Therefore, the operator on the left-hand side of cq. (15) is hermitian and so, the expectation value on the left hand 
side ought to be real. This is easily seen for the spin matrices on the right-hand-side as R&tyco/T) = cxp[a; • J/T] is 
hermitian. Showing that the integral on the right-hand-side of eq. (15) is real requires a more elaborate calculation 
which is carried out in Appendix B. 

Secondly, it should be pointed out that the matrix element in (9) has a spacial integral form which allows us to 
readily write a phase-space distribution: 

/(x,p) Tff = Ae- £ / T e ix -("- R ^»/ r )" 1 ('')) ^(^(W-'Ra^/TJWJ + ^CW^^/rJWt- 1 ))^ (17) 

where A = exp[^g/T] is the fugacity. Unlike its integral over the system's region, the function in (17) has, strictly 
speaking, no straightforward physical meaning and it does not need to be real or positive. Eventually, one could could 
calculate from it some more meaningful function like the Wigner's function. However, for very large systems, one 
expects the function (17) to become the actual phase-space single-particle distribution. For a cylindrically symmetric 
system to be macroscopic, the product TR, being R its maximal transverse radius, is supposed to be large, namely: 

Ti?> 1 

Hence, since a relativistic bound ujR < 1 is implied in eq. (2), for proper macroscopic systems (and indeed for most 
practical purposes) the ratio between lj and T is small, i.e.: 

& « • <»> 

where we have purposely restored the natural constants. Consequently, the difference between momenta in the 
exponent of (17) can be well approximated by the lowest order term in oj/T, i.e.: 

P — Ru)(*w/T) _1 (p) = p - cosh ^ p — i sinh ^ Cj x p + (1 - cosh ^) p • ww ~ v^Q> x p (19) 

so that the phase-space distribution function (17) reduces to: 

/(x,p) Tff = Acxp[- £ /T-x. ( W x p)/T]i (D s (\p}-^(iw/T)\p}) + D s (\pYRz(iw/T)\p}^)) Ta 

= A exp [s/T + p ■ ( W x x)/T] \ (D s ([p^R^ /T)[p\) + D s (^^(iu/T)^- 1 )) r<j 

= Aexp[- £ /T + p.v/T]i(i) s (b]- 1 R^(Wr)b]) + ^ S (b] t R^(Wr)b] t - 1 )) rff (20) 

where we have used eq. (2). The eq. (20) is the unnormalized single-particle phase-space distribution for the ideal 
rotating relativistic Boltzmann gas, seen by the inertial frame, and it turns out to be a non-diagonal matrix in the 
spin coordinates. The eq. (20) is one of the main results of this work and clearly shows that for a rotating gas the 
polarization states are not evenly populated. Furthermore, if the total number of particles is not fixed, (20) coincides 
with the density of particles in phase-space, including spin degrees of freedom. In order to obtain the properly called 
phase-space density in (x, p) we have to take the trace of the matrix (20): 

/(x,p) = £/(x,pU = Ae-^+P-/^) (21) 

a 

being: 



(22) 
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The exponent in (21) can be rewritten in a covariant form by introducing the temperature four-vector: 



R 1 ( \ 1 1 



(23) 



being 7 = (1 — v 2 ) -1 / 2 . Thereby, the invariant local temperature To = 1/ ' \[W measured in the comoving frame differs 
by a factor 7 from the global temperature measured by the inertial observer. Since 7 increases as a function of the 
radial distance, in a rotating gas at equilibrium there will be a steady gradient of local temperature; since also [ijT 
must be constant, there will also be a steady gradient of the local chemical potential, i.e. the local chemical potential 
will be given by 



Mo = 



(24) 



both these effects have been pointed out by Israel [13]. Finally, the phase-space distribution can be written in the 
compact form: 



/(x,p) rCT -Ac-^i(^(b]- 1 R^^/T)b]) + ^ s (b]tR (i (z W /T)b]t- 1 )) i 



(25) 



which definitely shows that the ideal rotating gas is made of cells moving at the four-velocity u{x) pertaining to a 
rigid motion. 

The rotational grand-canonical partition function of the ideal Boltzmann gas can be readily calculated by using the 
single-particle matrix element of eq. (15) for a = t, integrating over momenta, summing over polarization states and 
finally summing over all possible particle multiplicities N, dividing by N\, as usual. This eventually yields the simple 
result, in the approximation (18): 



exp 



^d 3 x yd 3 p/(x,p) =cx P A X (^)^d 3 x yd 3 p 



-/3(x)-p 



(26) 



III. THE IDEAL ROTATING GAS AS A PERFECT FLUID 



We have shown in the previous section that the globally equilibrated rotating gas can be seen as a set of locally 
equilibrated cells moving with a rigid velocity field (2), at least in the limit lo/T <C 1. In this section, we will reinforce 
this interpretation providing the relevant stress-energy tensor and the local thermodynamic relation. 

The stress-energy tensor T^ u can be calculated directly from particle phase space density /(x, p): 



d 3 P 



j/y7(x, P ) 



thus, for the ideal relativistic Boltzmann gas, according to (21), and using the four-temperature vector /3: 



t^ = xy(-')—— [ ^e-e-> 

X \T) dp dp J e 



(27) 



(28) 



The integral on the right-hand side of (28) is apparently Lorentz-invariant, therefore it can only depends on yT^ 2 
1/T ; if we define it as F(yJ~fl 2 ) we get, after some algebra: 



= A X (£) F" ( ^)^+Axg) nVP) 
= Axg) F'W)^ + Axg) ^gV 



g^P 2 - p»p v 



(29) 
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We are therefore led to identify the pressure and the proper energy density with 

p = _ Ax «M = T / d 3p /(x , p) 



2 



T 



p = A X (£)f"(v^) = I d 3 p( £ -vp)/(x,p 



) (30) 



and the stress-energy tensor thereby becomes: 

T» v = (p + p)uV -pg^ (31) 

which is the required form for a fluid at full thermodynamical equilibrium [13]. The identification of pressure with 
the expression in eq. (30) dictates that u/T = -yuj/T is a Lorentz invariant quantity; this will be further discussed in 
Sect. V. 

The correctness of the identification of pressure and proper energy density is confirmed by the analysis of the basic 
thermodynamical equation. For a general system in the rotational grand-canonical ensemble, entropy can be obtained 
from its definition: 

S = -tr[p w log/5^] 

By using the expression of p u in (7), we readily arrive at the basic thermodynamical relation for this ensemble involving 
the expectation values of global quantities: 

TS= (£?)- /x<Q)+T log Z„-w<J) (32) 

where T is the global temperature. The equation (32) can be transformed into a local relation by taking into account 
that \ogZu is in fact a spacial integral over the system's domain V. Comparing (26) with the first of (30) we can 
write: 



logZ. 



1 J d 3 xp(x) (33) 



This equation has been derived here for the ideal Boltzmann gas, but it holds for a general rotating system at 
thermodynamical equilibrium, as we will see. The mean values of energy, charge and angular momentum can be 
obtained from the partition function (8) by taking derivatives with respect to temperature, chemical potential and 
angular velocity. As the operators Q, H and J commute with each other, this gives: 



;H)-/i<Q) + w<j)) 



a log 5, = j_ 

dT T 2 
dlogZ u _ 1 ~ 

dp ~t w; 

d\ogZ u= 1 ~ 

doj T x 1 v 1 

Thus, if log is an integral over the system's domain of some function, the same happens for its derivatives and one 
can then re-write the entropy as an integral: 

r S = r/ d 3 x ^=/d3 x ^-,« + 15|^-„.f) (35) 

J v d J x J v d d x d J x d d x d J x 

The global temperature T is independent of x and one can then assume that the above relation holds locally: 

dS _d(E) d(Q) dlog^ d(J) 
T d5x-^r _M ^x" +T ^x __a; ' d^x" (36) 



2 The first equality in (30) can be easily proved by taking the derivative of F(y 
proved by taking into account that 



7J 2 ) with respect to 0°, while the second is most easily 



d 

rw m 1 ~ d^ 
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In general, the angular momentum density can be written as the sum of an orbital and a spin part: 

f>= x xf> + f> (37) 

d^x d^x d d x 

where d(P) /d 3 x is the linear momentum density. For the present, such decomposition can be explicitely demonstrated, 
for the ideal relativistic Boltzmann gas, only for the component of the angular momentum along u>. Indeed, plugging 
in the eq. (34) the partition function (26) and using the (23) and (2): 

^ogZ u _ /■ 3 l ^ tUr ^ , X '(t) 



(J • &) = T—^ = J d 3 x « • J d^p (x x p) /(x, p) + /(x, p) (38) 

which allows us to clearly identify the orbital and spin part of the angular momentum density. The complete calculation 
of all the components will be presented in Sect. V. 

By using the eq. (37) and recalling the relations T = T0/7, p = po/j, the (36) can be rewritten as: 

m 1 dS ld(E) . A d(P) ld(Q) ldlogZ w 1 d(S) 

7 d J x 7 d d x 7 d J x 7 d^x 7 d J x 7 d J x 

where we have multiplied and divided by 7 = (1 — v 2 )^ 1 where needed. Now: 

^~tt^ ~ 7( w x x ) • ^rr^ 

d d x d a x 

is manifestly the energy of the local cell with volume d 3 x seen in the cell's comoving frame, i.e. its mass density, so 
that we can also write the (39) as: 

IdS ld(M) ld(Q) l dlog^ ld(S) 

-M3 "To = To ^0--T^ ^^0 7 W T3— (40) 

7 d J x 7 d d x 7 d^x 7 d J x 7 d J x 

In the eq. (40) now appear densities rescaled by the Lorentz factor 7, i.e. densities measured in the comoving frame 
or proper densities. We can therefore write the above equality in the more familiar form: 

T s = p - Li q +p- juj ■ er (41) 

where s is the proper entropy density, p the proper energy density, q the proper charge density and a the proper spin 
density. The term Td log Z LU /d 3 x has been identified with the pressure on the basis of (33), which has been proved for 
the ideal Boltzmann gas. Nevertheless, since only through this identification the eq. (41) gets the form of the basic 
local thermodynamical relation, the (33) must have a general validity. 

With respect to the usually known form, the (41) contains an additional term owing to the non-vanishing spin 
density. In the non-relativistic limit this term becomes lj ■ a and is in fact responsible for the Barnctt effect [14], 
i.e. the spontaneous magnetization of a body spun around its axis when reaching thermodynamical equilibrium [15]. 
Unlike all other terms, the spin-rotation term has a non-covariant form. Expressing it, as necessary, in a fully covariant 
form, requires the introduction of two tensors, the acceleration tensor fi and the spin density tensor a. It will then 
become clear that the general covariant form of (41) is: 



T s = p-fi q+p- -fl^a^ (42) 



Deriving the tensors and a will be the subjects of the two forthcoming sections. 



IV. THE ACCELERATION TENSOR 



What is the right relativistic generalization of the us vector in a rigid velocity field (2)? This is the acceleration 
tensor SI, which is discussed for instance in ref. [16], as well as in several papers quoted in the bibliography, e.g. [5]. 
For the ease of reading, we spend this Section to introduce and describe the main features of this tensor. 

The starting point is the definition of an observer and a moving frame in the (Minkowski) spacetime, which is 
meant to be a generalization of an accelerated frame in non-relativistic physics. An observer is simply a world-line 
in spacetime. To this world-line one can attach a set of four unit vectors, called a tetrad, which are a basis of the 
tangent space in each point of the world-line are in fact a reference frame for the observer. The four tetrad vectors ej 
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are labelled by indices with the convention that the time-like vector eo is the four-velocity u, i.e. the vector tangent 
to the world-line. The other three space-like unit vectors can be chosen at will. The tetrad is orthonormal in the 
tangent space, thus: 

ei ■ ej — e i e^j — rjij 

being r\ — diag(l, — 1, — 1, — 1). Furthermore, any vector of the tangent space can be decomposed along the tetrad, 
hence the tetrad is complete: 

I = ^ ei <g> e i e in 4 = 0„„ (43) 

i i 

Since the tetrads at different space-time points of the world-line are orthonormal, the evolution of a tetrad along the 
observer world-line can be described by means of the Lorentz transformation A connecting the initial tetrad at the 
proper time To with that at the proper time r. Taking advantage of the group structure of Lorentz transformations, 
it is also possible to write the tetrad at the time r + dr as: 

e,(r + dr) = A(r + dr)(e,(0)) = A'(dr)A(r)(( ei (0)) = A'(dr)(( ei (r)) 

where A' is the infinitesimal Lorentz transformation connecting ej(r) to e^r + dr). Using the generators of Lorentz 
group J^: 



e;(r + dr) 



((e i (r))=e i (r)-ifi'"'drV(( ei (r)) (44) 



where, by definition fl^dr are the parameters of that infinitesimal Lorentz transformation. Reminding that: 

(V)£ = * (<£swj - Kg,*) 

and that W is an antisymmetric tensor, we can then rewrite (44) as: 

= ef = -^""(V)?ef = Sl af) e ip (45) 

or, in a more compact way: 

e 4 = • e, (46) 

Although we are working in a flat spacetime, we emphasized with symbols needed for general curved spacetimes, that 
the derivative in (45) is indeed the covariant derivative of the tetrad vector along the world-line, what is also denoted 
with an upper dot. If we are dealing with a fluid, then we will have to consider a tetrad field ej(x), because to each 
hydrodynamical cell will be associated a world-line. In hydrodynamics, the derivative along the cell trajectory is 
called the convective derivative and can be expressed as u ■ d. 

The tensor f2 in eq. (45) is the correct generalization of the angular velocity and acceleration and will be henceforth 
referred to as the acceleration tensor, following ref. [17], as has been mentioned. The eq. (45) generalizes the classical 
equations of motion of the axes of a rigid frame in space, i.e.: 

de ,„_. 
-=.xe (47) 

We will make the relation between and u) explicit for the uniform rotating fluid shortly; for the general relation, 
see Appendix C. The equation (45) can be easily inverted to obtain an explicit expression of W . For this purpose, 
one has to multiply both sides of (45) by e 4/3 , sum over the index i and take advantage of the completeness relation 
(43). The result is: 

W = (48) 

i 

In this form, the tensor ft is not manifestly antisymmetric, but this can be readily shown by using again (43): 
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We can now move on to determine the acceleration tensor for the rigidly rotating fluid. Yet, before doing this, we face 
another fundamental question, that is how to attach a tetrad to a hydrodynamical cell. When introducing the tetrad, 
we have invoked the world-line of an "observer" , which is always tacitly associated with an inner reference frame made 
of three orthogonal unit vectors, whereas for a material fluid pointlike cell, such association is not so straightforward 
and could be done in several manners. For instance, it can be done by using purely geometrical considerations, like 
e.g. Fermi- Walker transport or considering the instantaneous Lorentz boost without rotation connecting the tetrad to 
the inertial frame axes; or, rather, it can be specified by directions relevant to actual physical quantities like e.g. the 
spin density tensor discussed in the next section. In general, we would like to have the most natural definition of the 
tetrad, possibly determined only by the four-velocity u{x) without invoking additional physical fields. Such a natural 
tetrad exists and it is provided by the Frenet-Serret unit vectors associated to the world-lines of hydrodynamic cells 
[18], i.e. the four-velocity field tangent lines. Even if this definition depends solely on the field u(x), it is not unique 
because, at it is well-known, the four-velocity can be associated with the motion of charge (Eckart frame) or energy 
(Landau frame). Notwithstanding, in the present case of full thermodynamical equilibrium, all velocities must be 
the same [13] and for a macroscopic fluid with finite angular momentum are given by (2). We will later show, in 
Sect. V, that this tetrad is precisely the one which makes the spin density tensor, at least for an ideal Boltzmann gas, 
proportional to the acceleration tensor. Other choices would result in an acceleration tensor not proportional to the 
spin density tensor. 



CO 




FIG. 1: Picture of the rotating gas and the local Frenet-Serret reference frame or tetrad. 



For a uniformly rotating fluid, the Frenet-Serret vectors can be easily determined and are those corresponding to 
the axes of a rigidly rotating frame, that is, in the inertial frame components (see fig. 1): 

e = u = (7, 7V) = (7, jvj') ei = (0, i') e 2 = (jv, 7/) e 3 = (0, k') (49) 
whence we can calculate the derivatives along the world-lines. For i = 1,2, 3: 

^7^(0,7^ ei ) (50) 

because 7 and v are constant along the field lines. Therefore, the Frenet-Serret tetrad is such that the time derivative, 
in the fixed inertial frame, of the space vectors yields the simple law (47) involving the assigned u> vector. Taking the 
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derivative of eo, one gets the acceleration four- vector A: 

^ = ^^A = 1 ^ = ( 7 4 v • a, 7 4 v • a + 7 2 a) = (0, 7 2 a) = (0, 7 2 u, x v) (51) 

where, by definition, a = dv/di and we have used the explicit form of the velocity field (2) implying a = uj x v. We 
are now in a position to calculate the acceleration tensor from the formula (48). The easiest method is to multiply ft 
scalarly on both sides by the basis unit vectors of the inertial frame, that is: 

fV =E li -^2e i <g>e i -E v (52) 

i 

which leads us to, being uj = u)k = uik.', and using eqs. (51) and (50): 

/0 0\ 



ft 



jui 
- 7 w 
\o 0/ 



(53) 



It is worth pointing out that the double contraction of ft with itself yields: 

ft : ft = fl^W = 2 7 2 w 2 (54) 

and, therefore, 7 w must be a Lorentz invariant; this is further illustrated in Appendix C. If 7 w is a Lorentz invariant, 
so is the ratio ui/T as: 



^'■"-Jira-.a (55) 



ft ■. ft 



T 7 T T V 2 

As a consequence, pressure in eq. (30) can be written in a manifestly covariant form: 



p = A^^ X (y^ 2 ft:ftj (56) 
as well as all other local thermodynamical quantities. 

V. THE SPIN TENSOR 

We now come to the central subject of this work, i.e. the determination of the spin density tensor for the ideal 
rotating gas and its relation with the acceleration tensor. Some preliminary definitions are necessary. For a fluid, 
angular momentum conservation is expressed by an equation like: 

dxJ X fU/ = d x (S x ^ + x»T Xv - x»T x ^) = (57) 

where S X '^ V is a rank 3 tensor antisymmetric in the [iv indices commonly known as spin tensor. In classical field 
theories, the spin tensor can be made vanishing by a suitable redefinition of the stress-energy tensor, well known 
as Belinfante symmetrization procedure. However, such procedure does not leave the thermodynamic relation (42) 
invariant even when ft = 0, because both energy density and pressure are affected by the Belinfante transformation 
(see Appendix D for a detailed calculation) and the transformed entropy would have a different local and global 
value, which is unacceptable. The profound reason of the lack of the freedom of changing spin tensor and stress- 
energy tensor in presence of thermodynamics, is the fact that entropy is an effective counting of degrees of freedom, 
a concept tightly related to the discrete and quantum nature of matter and definitely outside a purely classical field 
mechanical framework. 

At full thermodynamical equilibrium we must have a symmetric stress-energy tensor in the form (31) and therefore 
the spin tensor must be conserved as well: 
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Furthermore, to meet the requirement that all currents must be parallel at full thermodynamical equilibrium, the 
spin tensor must be factorizable into an antisymmetric spin density tensor a and the four-velocity: 



(58) 



To determine the spin density tensor for a rotating ideal Boltzmann gas, it is necessary to calculate total angular 
momentum density and subtract the orbital part. A separation between angular and spin part has already been 
encountered in eq. (38) for the component of the total angular momentum along u>. The orbital angular momentum 
density must look like: 



dL 

d 3 



[ r 

-= / d 3 p (x x p)/(x, p) 



hence, if such a form is retrieved in the calculation of the total angular momentum, the rest can be unambiguously 
identified as the spin density. In the particular case of the Boltzmann ideal gas, particles can be handled independently 
of each other. This makes it easy to calculate the spin density tensor, what can be done by simply summing up single 
particle contributions. We then need to calculate the expectation value of the angular momentum tensor (i.e. the 
Lorentz group generators) over single particle states, that is: 



d 3 p (p, a\f v exp[(-/i + nq + u ■ j)/T]P v \p, a) 



(59) 



where the normalizing factor 



(N)=J^d 3 x J d 3 p/(x,p) 



concides with the mean number of particles (see Sect. II). We will proceed separately for the space components of the 
angular momentum tensor, i.e. the properly called angular momentum J, and the time components, i.e. the boost 
generators K* = J°\ Taking the z axis as the w direction, we can write the mean values in (59) as derivatives: 



«J*» 



— y 



d ~ - - 

d 3 P i^riP, °1 R ei i<Pi) exp[(-/i + fiq + ujj z )/T]P v \p, a) 
o<Pi 



01=0 



((%)) = J^Y.J^V^iP^lU^cMi-h + ^q + ^^mPvlp^} 



6=0 



(60) 



where Rg, (</>;) = cxp[— i<pili\ is a rotation of axis e; and angle (pi and Lg ; (6) = exp[— i£;K/] is a Lorentz boost of 
direction e; and hyperbolic angle The matrix elements in the above equation can be calculated in a similar fashion 
as for that in eq. (9). One just needs to insert a resolution of the identity made of momentum-spin eigenstates between 
the operators Rg ; (</>;), 1-^(6) and the density operator and work out the obtained expression following the procedure 
drawn in Sect. II from eq. (9) to (15). This leads to: 

((ji)) = ^E/ d3 P c ~ £/TWT ^/ dSxe^P-^^)" 1 ^)- 1 ^)) 

x \ (^ s (W- 1 Re,(0OR^(Wr)b]) + ^ s (M t Re i (0OR^(Wr)b] t - 1 )) (T(T a 

= _L / d 3 p e -e/T +m ,T i 9 f d 3 xe4 x.(p-R li ^/T)-iR, ! (0 ! )- 1 (p)) trjD S (R . i( ^ )R(i)( , cj/T)) 

(A) J dcpi J v 



0!=O 



(61) 



and 



((h)) 



= % 



d_ l_ 

96 (AO 



d 3 x e^Ui^iv-^^/T)- 1 ^ {ii)- 1 { P )Y 



x i(^([ P ]- 1 Le ! (6)R^^/r)[ P ]) + ^(WtLe ! (6) t - 1 R^^/r)b]t- 1 )) ( 



6=0 



" Z (N)S 



d 3 pi 



-U, (iirHpf/T+M/T 



1 (p)°Sv & 



x e 



x 1 [tr^ s (Le ! (6)R^^/T))+tr J D s (Le ! te) H1 R^(^/7 1 ))] 



6=0 



(62) 
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Note the difference between the traces of the finite dimensional representation matrices in (61) and (62) which is 
owing to the non-unitarity of the SL(2,C) representation of boosts, i.e. D s (L^))^ 1 ^ D S (L(£)). 

Let us first focus on the eq. (61). If we let p' = R 6l (<fii) _1 (p), the energy as well as the momentum integral measure 
does not change and the derivative with respect to <f>i can be worked out: 



(0j» 



,\- [ d 3 p' C- £ '/T+ M /T z 9 f d 3 xc4 x.(R, ; (0 ! )(p'-R^(-/T)- 1 ( P ')) tr ^S (R6!((/);)R( . ( ^ /r)) 

(A) J o<pi Jv 



1=0 



J d 3 p' e - £ '/^/T J d 3 x e iMp')-^/T)- l ( P ')) [ (x x p /) . e lX (^) + tr£ s (J ; R^/T))] (63) 



(N) J jv 

where we used the (22). The (63) can be rewritten in a more compact form recalling the calculations of Sect. 2: 
( ^ )) = MX d3x / d3 P/( x 'P) (xxp) + -^ytr[£> s (JR«i(iw/r))] 



(64) 



where we have redefined p' as p. In this form, eq. (64) allows us to pin down the spacial part of the spin density 
tensor which is the second term in the integrand of (64) and more specifically: 



dS 

d% 



= / d 3 p f( x ,p)-l—trD s (JR d ,(iLu/T)) 



The z component of the trace is readily calculated as: 



tvD s (i z R^/T)) = —tvD s (R^/T)) = ^tr exp[^ 5 (J z )/T] = X 'g) 



du/T 



(65) 



(66) 



As far as the x, y components are concerned, we first note that for any analytic function /(J) of the rotation generators, 
or a representation thereof, and for any operator such that [0, J z ] = 0, we have: 



tr(J Xil/ /(J z )0) =0 



(67) 



To prove it, we just calculate the trace of S XtV S^ +1 0: taking advantage of the ciclicity of trace and the Lie algebra 
properties: 



tr (J^jf +1 0) = tr (j XiV J^OJ z ) = tr (J J^jf 0) 

= tr (J^Jjf 0) ± itr (J^jf 0) = tr (j^jf +1 0) ± itr (J^jf 0) 



(68) 



Comparing the first and last side of the above equation, we conclude that tr[J yj2 J^ 0] = for any N and therefore 
the (67). We now apply the above reasoning to our specific case where /(Jz) = cxp[cjJ z /T] and = 1, which both 
fulfill the aforementioned conditions and thus conclude that: 



\xD s (l XtV Ra(iu/T)) = 



(69) 



As a consequence, the spin density (65) turns out to be, in view of (66) and (69): 

^ = y dp/(x ' p) ^y w=n(x) ^y w 

being n(x) = J d 3 p /(x, p) the particle density. Similarly, the proper spin density (see Section III) turns out to be: 



(70) 



' ■ '■ ■-■ • , \X ( t) ' "X- ( T ) - 

er = — — — = —nix) — j^-i-OJ = v — t+t-U) 



1 dS _ I 

7 d 3 x 7 



being v the proper particle density: 



I I f T 

v=-n{x) = - / d 3 p/(x,p) 

7 7 J 



r J d 3 p /(x, 



P) 



T 



(71) 



(72) 



where we have used the equation (30). We note in passing that the x, y components of the orbital part of the angular 
momentum density vanish as well because the system (and the phase space distribution) is geometrically invariant by 
rotation around the z axis. 
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We can now turn to the eq. (62). We have seen, for angular momentum, that the derivative with respect to group 
parameter produces two terms: the first, obtained by taking the derivative of the spacial integrand, gives rise to an 
"orbital" part, while the second term stems from the derivative of the matrices of the finite-dimensional SL(2,C) 
representation and gives rise to the actual "spin" part. Similarly, for boosts, the derivative of SL(2,C) representation 
matrices yields the contribution to the spin tensor density while the derivative of the spacial integral generates an 
orbital part: 



«**» 



orb 



(TV) 



UtfO-W Jv 



(n) Xx (t) 



U.tfO-W Jv 



d 3 x e^^i^b-MWr)- 1 !-^,)- 1 ^)) 
d 3 x e i £?=i x'-(p-R tf( (ia,/r)- 1 L. I K,)- 1 ^)) 



6=0 



To go on, we need to work out the action of a Lorentz boost of hyperbolic angle £ along the direction n on the 
four-momentum p: 



Lfi(0(p)° = P°cosh£ + n-psinh£ 

U((,)(py =P { +p°smh£h l + (cosh£- l)n-pn l 

Let p' = Le,(£;) -1 (p) and implement this change of variable in the integral of eq. (73): 



orb 



)-R li (^/T)- 1 (p'))* 



6=0 



(74) 



(75) 



The derivative can be worked out by using (74): 



((h)) 



orb 



/3/2 



E '/ T / d 3 xe i 2=' / ( u .^'< !, ''- R »(™'' T '" 1 < ,, '»' 



6=o 



3i , 



2<r 



7 p ■ e; - £ X • e; 



^.(p'-R^/T)- 1 ^')) 



W) I d 3 X / d^p 



3i „ 

— p • e ; - ex • ei 

2e 



/(x,p) 



(76) 



where, in the last expression, we have redefined p' as p. It is now easy to realize that, in a rotating system, the last 
integral vanishes. In fact, we have: 



d«k» 



orb 



d 3 P 



= M§ p // x/(x ' p) -M £ X d3xx/(x ' p 



(77) 



The rightmost integral yields the origin of the inertial frame, which is centered on the axis of the rotating system 
according to the (2). The leftmost integral, on the other hand, is imaginary and the integrand is proportional to 
particle velocity; once integrated over momentum, this term yields the mean particle velocity which vanishes. Thus, 
for any constant C: 



k>U = j^J v d 3 x J d 3 p (CH _ ex)/(X; p) = o 



« 



On the other hand, the "spin" part of the eq. (62) can be written as: 



((h)) 



spin 



(N)S 
d 1 



d d pe~ 



7 d3 

Jv 



, i,, „-e/T+iiq/T / H 3 X ^-(p-R* (tw/T) 1 (p)) 



x ^ (tr [D^l^R^ioj/T))] + tr [D s (L^-^iuj/T))]) 



6=0 



(KiR&(iuj/T)) + D s {K\R*(iuj/T)) 



- w>X d3x / dSp/(x - p, 5 tr 

By using the known properties of representations, the last matrix can be written as: 

£ S (K ; )MWT)) + D s {K\R*{iu/T)) = D s {Ki)D s (R £t {iu/T)) + D s '(K})D S '(Ra(tw/T)) 



(78) 
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but in fact, for finite-dimensional SL(2,C) representations D S (K) is anti-hermitian, hence D S (K^) = — D S (K) and the 
above expression vanishes. Consequently, the contribution of boosts to spin density in a rotating ideal Boltzmann gas 
is zero, i.e.: 



d5 w 



= 



dK s 



or 



= 



d 3 x d 3 x 
Collecting the eqs. (71) and (79) we can finally write the spin density tensor a as: 



(79) 











o\ 








v m 








x(t) 








\o 











(80) 



The result (80) has three major implications. Firstly, for an ideal relativistic Boltzmann gas at thermodynamical 
equilibrium, the spin density tensor (80) is proportional to the acceleration tensor (53) defined with the Frenet-Serret 
tetrad, through a Lorentz scalar t depending on the intensive invariant parameters T , no and ^ : 



with: 



i = v 



1 x'(f) 



(81) 



(82) 



where, in (82), we have used the eqs. (54), (55), (72) and (56). The last form of i emphasizes its nature of Lorentz 
scalar depending on the Lorentz scalars /zq, /3 2 and f2 : f2. Also, since co C T and using the (22): 



v S(S + l)u v 5(5 + 1) v 5(5+1) 



70; 



IT 



2 v_ 5(5 + 1) 

KTn 3 



(83) 



where in the last expression we have restored the natural constants. 
Secondly, the basic thermodynamical relation (41) can be written as: 

1 



Tqs = p - poq + p- -\l:a = p- p, a q + p 
as anticipated in eq. (42). For an ideal Boltzmann gas the additional term is: 

uj 2 5(5 + 1) 



2 iL iiv<J 



2 



-iQ. : n 
2 



2 x(t) 



h'uj 2 5(5 + 1) 
~KT 3 



(84) 



(85) 



using (82). Once again, in the eq. (85), we have restored the natural units, which allows us to conclude that 
for ordinary macroscopic gases the correction to thermodynamic relation owing to spin is extremely small, being 
Tiuj/KT <C 1 under all circumstances. For instance, for a mole of a monoatomic gas with 5 = 1 at 273 °K, p = 10 5 
Pa and uj = 100 Hz, fiuo/KT = O(10" 12 ) and the energy density associated to spin-angular velocity term is of the 
order of 10~ 19 J/m 3 . However, in extreme situation, like e.g. in astrophysics or relativistic heavy ion collisions, this 
term could be much more important. Also, as has been mentioned at the end of Sect. Ill, the term coupling spin 
and angular velocity-acceleration is responsible for the magnetization observed in rotating bodies, the Barnett effect 
[14, 15]. 

Finally, the projection of the spin density tensor over the four- velocity field, is non- vanishing and proportional to 
the four-acceleration A through the factor i: 



Particularly, in the inertial frame, this vector has time and space components: 



(86) 



t = (0, rfu) x v) 



(87) 
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(see eq. (51) and it is thus orthogonal to both u; and the velocity. This result is in contrast with common assumptions 
in the hydrodynamics of fluids with spin, as has been mentioned in the Introduction, but it is not surprising as there 
is no principle reason for it to vanish in a system of many particles [3]. Indeed, while the vector t for a quantum 
relativistic particle in a pure state vanishes, this is generally neither the case for a many-particle states nor for a single 
particle in a mixture of states (see Appendix E). 

We end this Section by noting that the calculations herein performed provide an alternative way of determining the 
polarization four- vector of particles, which was computed explicitely for different spins in ref. [8]. The relevant proof 
is shown in Appendix F. 



VI. VIEWPOINT OF THE ACCELERATED OBSERVER 



Since the ideal gas we are dealing with is a macroscopic system, we can ideally attach to each point a physical 
observer (i.e. a tetrad) moving along the velocity flow lines. For him to observe local thermodynamical equilibrium, it 
is necessary that all (local) physical quantities are constant in his frame. If this is the case, and assuming the locality 
hypothesis [17, 19], stating that the accelerated observer measures the same quantities as the comoving observer 
(i.e. the inertial observer instantaneously coinciding with the accelerated one), the relation he finds among the local 
thermodynamical quantities are actual equations of state. Particularly, the relation (81) linking a kinematical quantity 
like the acceleration tensor and the spin density tensor is to be promoted to the rank of equation of state for a system 
in local thermodynamical equilibrium. 

The condition of steadiness for the non-inertial observer means, for a scalar: 

where t is the proper time. For vector quantities, this turns into a condition on the components with respect to the 
observer's tetrad: 

<LY* _ djX ■ e Q _ Q 
dr dr 

implying that the four- vector X varies according to: 

x = J2( x ' te >) = Yl x * ei + Y x * ii = Y x ' ki = Y xtfl - e * = n-x 



Similarly, for rank 2 tensors, one should have: 



dT'J 





dr 
implying: 

t = ^(T^e- ® ej ) = T3 (e« ® e-i + ^ <g> e 3 ) = ^ T ^ ( fi ' e i ® ej + e t ® ft ■ e,) = ft • T + T ■ fl T 

ij ij ij 

where T stands for the transpose. Applying the above equation to the acceleration tensor itself, one obtains that 
the condition of steadiness amounts to enforcing that the convective derivative of ft at thermodynamical equilibrium 
must vanish, i.e.: 

ft = (88) 

For a rigidly rotating fluid, all above conditions are fulfilled provided that the tetrad is the Frenet-Serret one. The 
local temperature T , the local chemical potential /i and the acceleration tensor have vanishing convective derivatives, 
as it turns out from eqs. (23), (24) and (53). The spin tensor density is also stationary in the non-inertial frame because, 
for the ideal gas, it is proportional to the acceleration tensor. We stress that this does not happen for other tetrads, e.g. 
a Fermi- Walker transported one; in that frame, the acceleration tensor components and, consequently, the spin density 
tensor components would vary and therefore, one could not properly talk about local thermodynamical equilibrium. 
This is why, in a rigidly rotating fluid, the Frenet-Serret frame is a privileged reference. 
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Comparing the view of the inertial observer with that of the non-inertial one allows to identify inertial forces and 
helps elucidating the role of spin. To start with, it is useful to take the non-relativistic limit of phase-space distribution 
(20) seen by the inertial observer. At first order in v, for a gas of spinless particles, this reads: 



/(x,p) oc cxp 



P 



v p 



2mT 



T 



exp 



P 



2mT 



+ 



(w x x) • p 



exp 



P 



2mT 



+ 



uj ■ L 



(89) 



which, once integrated, yields the invariant partition function. At each point x the inertial observer has to integrate in 
momentum the integral (89) because he deals with a rotating gas and, in order to take angular momentum conservation 
into account, he has to calculate a rotational grand-canonical partition function with anu-L term [9]. Conversely, for 
the non-inertial local observer, in the non-relativistic usual view, the particles of the gas as living in external fields, 
the centrifugal and Coriolis fields. He starts writing a lagrangian: 



L = imv' 2 + ^m(u/ x x) 2 - mx • (u/ x v') 
with centrifugal and Coriolis potential energies. From (90) a hamiltonian can be derived: 



11 p 2 

H = -mv' 2 - -m(u/ x x) 2 = ^ mp • (u/ x x) 

2 2 y ' 2m v y ' 



where, from (90): 



p = mv — mu x x 



(90) 



(91) 



(92) 



Since in the non-relativistic limit u> — u>' , the non-inertial observer then writes a single-particle partition function in 
the point x: 



J d 3 p e~ H ' T = J 



p exp 



P 



2mT 



v P 

T 



(93) 



with the same temperature T, as there is no difference in the non-relativistic limit. The result is the same the inertial 
observer would get, according to (89) and this means that pressure, entropy density and all local thermodynamical 
quantities are the same for the two observers. Yet, the physical interpretation of the non-inertial observer is different: 
for his point of view, the coupling between uj and x x p arises from the presence of external fields, i.e. inertial forces, 
in his frame the angular momentum of the system is simply vanishing. 

It is worth pointing out that, from a relativistic point of view, the effect of centrifugal potential in the hamiltonian 
(91) is already contained in the dependence of temperature on the distance from the rotation axis. The rotating 
observers at different distances r measure different local temperatures To(r) and write the single-particle partition 
function of an ideal gas without external fields in their comoving frame, where /3 = (1/T , 0): 



J d 3 p e-"* = J d 3 p e - £ / T °W = J d 3 p e ^hir)T ^ m J d 3„/ 



exp 



e 
T 



ev 
2T 



The last term yields, in the non-relativistic limit and disregarding the mass term: 

y' 2 1 (w x x) 2 " 



m 



/dV 



cxp 



2T 



—m- 

2 



(94) 



(95) 



which eventually, after the change of variable (92), turns out to be equal to the right hand side of eq. (93). This result 
was as a consequence, from the point of view of the non-inertial observer, of the presence of an apparent gravitational 
field introducing a position-dependent temperature [20], which is still an inertial effect. 

Similarly, the non-vanishing polarization of particles that he measures will be attributed to an inertial effect, a 
coupling between the acceleration tensor and spin of the form: 



which is, in his view, an additional term of the hamiltonian to be taken into account when calculating thermodynamical 
quantities. Clarifying the nature of this energy, its relation with Thomas precession and spin-rotation coupling [21] 
and calculating thermodynamical quantities in stationary accelerated frames will be the subjects of another work. 
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VII. SUMMARY AND OUTLOOK 

We have studied the physics of the ideal relativistic rotating gas at full thermodynamical equilibrium, in the 
Boltzmann approximation. We have shown that, provided that it is sufficiently large, the gas is a rigidly rotating fluid 
with a non- vanishing spin density. We have obtained the expressions of the phase space distribution (see Sect. II): 

/(x,p) rCT = Ae-^i(^(b]- 1 R li (^)b]) + ^ S (b] t R^(WT)b] t - 1 )) r(T 

as a matrix in the polarization states. From this distribution, we have inferred the spin density tensor and found 
that, at equilibrium, is proportional to the acceleration tensor constructed with Frenet-Serret tetrad, i.e.: 

where t is a Lorentz-scalar coefficient depending on the other thermodynamical variables. Furthermore, we have 
recovered the generalized fundamental thermodynamical relation 

T s = p - noq + P - ^ CVct^ 

including an additional term involving the spin density tensor, already used in literature in fact. We have shown (see 
Appendix E) that this term cannot be made vanishing by using Bclinfante transformations of the stress-energy and 
spin tensors because this would imply a change in the total entropy and so the spin density tensor has physical reality. 

The condition of full equilibrium makes it possible to study a relativistic fluid without the prerequirement of a 
theory of perfect fluids with spin and therefore, it allows to carry out an independent test of assumptions made in 
the past about the dynamics of fluids with spin. Particularly, we have found that the projection of the spin density 
tensor onto the four- velocity 

is not vanishing, unlike commonly assumed [6, 7], in agreement with general arguments of special relativity [3] and 
envisaged by Bohm and Vigier [4]. Therefore, an observer comoving with the spin density needs two non- vanishing 
space-like vectors, not just one, to describe it. 

The non-vanishing of the four-vector t has a considerable impact on the general theory of a relativistic perfect 
fluid with spin, which needs to be revised from the formulation of Halbwachs [5] and Ray [22]. This revision may 
open a window on new and unexpected features of relativistic hydrodynamics and its phenomenological consequences 
in various fields where relativistic hydrodynamics is a fundamental tool, from relativistic heavy ion collisions to 
cosmology. 
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APPENDIX A - Projection onto localized states 

The meaning of the projector Py is easicrly understood in a first-quantization framework for a single spinlcss 
particle, where it can be written a sum over a complete set of eigenfunctions vanishing outside the region V: 

p v = j2\ k )(k\ 
k 

with: 

u fe (x) = (x,er|k,er) 

The functions Uk are complete in the sense that: 

J2 4(x'H(x) - <5 3 (x - x')Mx)Mx') 

k 

where 9y (x) = if x ^ V. Therefore, one has: 

(p'|Pv|p) = J>'I*X*IP> = d'Vc-^- x X Ufe (x') ? 4(x)e J P- x = |d 3 xe l ("'»' x (96) 

which accounts for the integral part of eq. (12). 

In a quantum-field theoretical framework, the projector Py can still be defined and its matrix elements calculated 
[10, 11] giving rise to the eq. (12) for single-particle states. Although the derivation of (12) in refs. [10, 11] was partly 
based on an unproved conjecture, the non-trivial spin structure of eq. (12) can be indeed justified by the requirement 
that the projector onto localized states commute with rotation operators: 

for any angle <fi. Therefore, by applying single-particle relativistic states on both left and right hand side: 

(pV|PvR ft (0)b,<7> =^(p',a / |Py|R fl (<A)(p),r)£> s ([R fi (0)(p)]- 1 R fi (^)[p])ra (97) 
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whereas: 

<pV|R a (#Py|p,a) = (p,a\PlR A (^\p',aT = ^(p, alPylRft^)- 1 ^), r)^^^^)- 1 ^')]- 1 ^^)- 1 ^])^, 

(98) 

The (97) and (98) ought to be equal, so if we define: 

M(p,q) a ^ = (p,a\P v \q,0) 

from (97) and (98) ensues: 

M(p', R ft (0)(p)) J D s ([R A (^)(p)]- 1 R ft (0)[p]) - R fi (0)[R A (0)- V^Wfe W)~V)) (99) 

The last equation tells us that the matrix M cannot be proportional to the identity, as far as the spin is concerned. 
In fact, the eq. (99) can be fulfilled if, according to eq. (12) : 

M(p,q) = F(p -q) 1 - [D 3 ^}- 1 ^}) + £> S ([p]t[ g ]t-l)] 

being F a scalar function. Thereby, the eq. (99) becomes: 

^(MtRnWiRnW-W -1 )^ ^-RnW-V))* 

(100) 

Taking advantage of the unitarity of Wigner rotations (see eq. (14)), both left and right hand side of the above 
equation yield the same matrix multiplied by two different scalar functions: 

f( p ' R fi (0)(p))^ (^(bT'RftWb]) + ^(bfRftWO^- 1 )) 

= \ (fl 5 (b']- 1 R 4 Wb]) + fl s ([p'] t RiiWb'] , " 1 ))f(p-R*W 1 (p')r (ioi) 

so that we are left with the simple requirement: 

F(p' - R„(0)(p)) - F(p - Rft^J-V))* 
which is fulfilled if, again according to eq. (12): 

*V - RfiW(p)) <* / d 3 x exp[zx • (p' - R fi (0)(p))] 
and if the region V is invariant by rotation around the axis n. 

APPENDIX B - Reality of phase space distribution 

We want to show that the integral: 

f d 3 x cxp [ix • (p - R^(io;/r)- 1 (p))] (102) 
Jv 

is real. First, we note that: 

p-R^uVT)- 1 ^) 

has vanishing component along uj. Therefore, taking oj as the unit vector of the z axis, we can write the above integral 
as: 

/ d 3 x exp [ix T ■ (p T - R^(iw/T)- 1 (p T ))] 
Jv 
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and, by using the explicit expression in (19): 



/ d 3 x exp ixr ■ ( (1 — cosh ^)pr ~ i sinh ^lixpT 
Jv L V 1 1 



The imaginary part of the last integral reads: 



d 3 x sin 



to , 



(1 - cosh— )x T • p T 



exp 



sinh — x T • (u> x p T ) 



(103) 



In order to show that this vanishes, we need to find a rotation of xt around uj, i.e. in the transverse plane, which 
changes the sign of x-r • px and, at the same time, leaves xr • {ui x px) unchanged. This transformation in the 
transverse plane exists because pr and (cl> x pt) belong to that plane and are orthogonal. Once this transformation 
is performed, and the relevant change of variable x — > x' in the integral applied, the integral (103) is tranformed into 
its opposite as the region V is symmetric by rotation around uj, hence (103) vanishes. 



APPENDIX C - The angular velocity vector in relativity 

To better understand the physical content of the antisymmetric tensor ft, it is useful to decompose it into two 
four-vectors through the definitions: 



1 

^ _ ~2 



-e^ p M vp u a (104) 



which can be inverted to give [16]: 



&nv = e f _ ll , pa cu p u cr + A^Uv - A u Ufj_ (105) 



Both four-vectors A and ui are orthogonal to the four-velocity and are therefore space-like. It is easy to verify, by 
using (48) and reminding that eo = u, that A is just the four-acceleration vector du/dr. In order to grasp the physical 
meaning of the four-vector ui, it is necessary to first look at motion of the tetrad when ui — 0. In this case, the 
equation of motion of the tetrad four-vectors read, from (45) and (105): 

ef = A p e 2 ■ u - u p A ■ e, (106) 

This is the well known Fermi- Walker transport of the vector ej. Now, with respect to the Fermi- Walker transported 
tetrad, any other comoving spacial frame will be seen by that observer either as fixed with respect to the Fermi- Walker 
one or "rotating" with an equation like (47) with the time t replaced by the observer's proper time r, that is: 

e' t = J x e' (107) 

where the prime denotes that these are spacial vectors seen by the moving observer and not by the inertial one. In 
the Fermi- Walker transported frame the four- vector co reads: 

1 1 3 1 1 3 

,,, — f Q u Pn a — f V* p iv p p n' J — e A v ii p ii a e 'S^p iv p p n' T 

^ 2 t J,u P a / * i — 2 t J,u P a 2 jUy ' 9<7 / ; i 

i=0 i=l 



1 3 

= -gWoE^e? (108) 

i=\ 

being u a = 8% in any comoving frame. The last equality, along with (107) implies that in the comoving Fermi- Walker 
transported frame 

LO = (0,U/) 

which makes it clear the meaning of the four- vector lo. Now, to get it expressed in the inertial frame, it suffices to apply 
a pure Lorentz boost (without rotations) to the components seen by the comoving observer, with an instantaneuous 
velocity v with respect to the inertial observer. Therefore: 

u° = 7cj'° + 7V • u>' = 7V • u>' 

uj' = lv h + + ^_J- W W = lj h + 1_1 V VV (109) 
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where we also used u/° = 0. 

The last step is to provide a definition of angular velocity vector us in the inertial frame. We could, for instance, 
define it as the spacial part of the lo four- vector in (109), but it would not coincide with the angular velocity vector 
to that has been defined for a uniformly spinning system related to the velocity with eq. (2). Instead, we will use the 
latter as definition of uj and thereafter derive its relation between u/. For a uniformly rotating system, we know the 
f2 tensor from (53) and this allows to readily determine the ui four-vector from (104): 

uj = (~/ 2 to • v, 7 2 u>) (HO) 

We will keep this as a definition of ui in general. One one hand, the (110), along with the (51), allows us to write the 
tensor in general as a function of the acceleration a and the angular velocity u> seen by the inertial frame, by using 
the (105). This leads to, in the inertial frame: 






7 3 a — 7 3 (w x v) 


— 7 3 a — 7 3 (w x v) 


7 3 u; — 7 3 u; • vv + 7 3 a x v 



(111) 



Finally, comparing (110) with (109) we obtain the transformation rule of the angular velocity vector between the 
inertial and comoving frame: 



1 7-1 1 1 
u> = — o w' + -^-5-5- v • u/v = —^oj' H -vw'v (112) 

7 Z j z v 7 7+1 



APPENDIX D - Non-invariance of thermodynamics under Belinfante transformation 

The conservation of energy-momentum and angular momentum for a continuous system can be written as system 
of two partial differential equations: 

d^ v = 

T being the stress-energy tensor and S the spin tensor. This set of two equations is invariant in form under a 
transformation [5]: 

where <I> is an arbitrary tensor of rank three which is antisymmetric in the last two indices. 
Choosing ^ x ^ v = S X '^ , the new spin tensor S' vanishes and the new stress-energy tensor 

T ^ v = T ^ + 1^ ^ s \ tflv _ _ S »M} ( 115 ) 

becomes symmetric; this is the well-known Belinfante symmetrization procedure. However, this transformation does 
not leave, for a system in full thermodynamical equilibrium, the basic relation (42) invariant, whether the : a/2 
term is included or not, and, moreover, it also changes the total entropy. To show this, we first remind the definitions 
of proper energy density and pressure for a given stress-energy tensor: 

p = T^u ll u v P=\{u ll u v -g IJlv )T» v (116) 

For a system at full thermodynamical equilibrium, such as a rigidly rotating fluid, the stress-energy tensor must have 
[13] the form: 

T» v = {p + p)u»u v - pg» v (117) 
as well as a spin tensor which must factorize into a spin density second-order tensor and the four-velocity: 
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The stress-energy tensor in (117) is already symmetric, yet the spin tensor is not vanishing. We could then think of 
applying the Bclinfantc symmctrization procedure described above in order to make the new spin tensor S' vanishing, 
thereby obtaining a new symmetric stress-energy tensor T 1 . Using eq. (115) and (118), this new tensor reads: 

T'^ = T» v + \d x (a^u x - o Xv u» - <r A V) = - \d x {<r Xv u») - \d x {<r x »u v ) + \a* w + \o» w d ■ u (119) 

Zi Z Z Z Z 

At thermodynamical equilibrium the last two terms vanish: the second one, because the velocity field is rigid, and 
the first one because of the (88) and the vanishing of convective derivatives of T , po- 

& = (til) = L(T , Mo, ^ : ^ + id = 

Thus, the (119) becomes: 

T'^ = - id A (cr A V) - ^dxia^u") (120) 

For the stress-energy tensor according (120), the new energy density and pressure read, applying the (116): 
p' = T'^u^Uv =p- Ufl u v dx{a Xv u^) =p- u v 3xo Xv 

P' = \{P -TT) = \{P- ^dxa x » - T; + dxi^uv)) = p + rf Xv dxu v (121) 
The thermodynamical relation (42) in the new variables should read, as S' — and assuming the invariance of q: 

T s' =p'+p' - p q' = p + p-fi q- u v d x o Xv + ^a Xv 8 x u v (122) 

Hence, in order to keep the thermodynamic relation invariant, the entropy should be modified by a term depending 
on the spin density tensor: 

T Q s' = T s - u v d x o Xv + \c7 Xv d x u v + Ulx v <J Xv = T s - d x {<y Xv u v ) + a Xv d x u v + ^cr Xv d x u v + Ul x „cr Xv (123) 

The last expression can also be rewritten as, by using the definition (86) and taking advantage of the antisymmetry 
of a: 

T s' = T s-8-t+ ^a Xv {d x u v - d v ux) + hlxu<r Xu (124) 

It can be shown that, for a rotating gas, all terms added to entropy density on the right-hand-side of the above 
equation are strictly positive. Using eqs. (87), (82) as well as (54) and (53), the eq. (124) can be rewritten as: 

2 

T (s' - s) = -V • (t7 2 u; x v) + -i^u) • V x 7V + lj 2 w 2 

o 

The vector field lj 2 uj x v is proportional to the centripetal acceleration and therefore has a negative divergence; the 
curl of the vector field V x 7V for the velocity (2) yields 2-fuj. Hence, all three terms on the right-hand are positive 
and total entropy should be changed by a finite positive amount. 

The latter condition is clearly physically unacceptable. Since the rigidly rotating gas is a system in full thermo- 
dynamical equilibrium, the new entropy density obtained after the Bclinfantc transformation should at least have a 
vanishing integral over the system's region, which is not the case. 

Therefore, we conclude that making the spin tensor vanishing is impossible and it must have a physical reality. 



APPENDIX E - The four-vector t 

The four-vector t, for a single massive particle, is defined as the expectation value of the operator t m = J^P^ 
divided by the mass. If the particle is in a pure state \p, <r) then the expectation value: 

f = -(p,a\r^\p,a) = -(p,a\J^P v \p,a) 
m m 
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transforms as a vector because so does the relevant operator. In other words, by using the very definition of single 
particle states from the theory of Poincarc group representations [12, 23]: 

f = l<pVlM _1 > v ilMbV) = ^(pVlM^^bV) = M^( P V|J A ^bV) (125) 

where p° = (to, 0, 0, 0). Therefore, to prove that t = we just need to show it for a particle at rest. In this case: 

(p ,a\j^P„\ P ,*)=m(p°,CT\J» \ P ,a) 

The boost generators however anticommute with the parity operator and since the states of a single particle at 
rest have a definite parity (which is, by definition, the intrinsic parity of particles) the right hand side of the above 
equation ought to vanish. Thence, in view of (125): 

t" = -(p,an,a) = -{p,a\J^P v \p,a)=Q (126) 

TO TO 

The validity of eq. (126) is readily extended to any pure state. First we note that the above reasoning can be 
extended to a matrix element with different spin components and same momentum: 

(p,<T\T»\p,a')=0 

For a general pure state, it suffices to decompose single particle states onto momentum-spin eigenstates and show 
that: 

(p',a'\T»\p,a)=0 (127) 
for p^p' . By using standard Poincare group commutation relations: 

[J A<T , P v \ = iP x g 1 " 7 - iP a g vX 

it is easy to show that: 

[f, P v \ = iP^P" - ig^P 2 (128) 

and therefore, if p ^ p': 

(p',a'\T»P»\p,a) = (p',a'\P^\p,a)=p v (p',a'\^\p,a) = p lv {p' ,a'\t»\p,a) 

that is: 

if -p' v W^\^\p,cj)=Q 

Since the last equation applies to all p's components, and p ^ p' , it turns out that the only possibility is that the 
matrix element of t vanishes, i.e. the eq. (127). 

Altogether, all single-particle matrix elements of the operator f vanish, i.e.: 

(p',a'\T»\p,a) =0 

for any p, p' , a and a', and thus the expectation value of f, hence t, vanishes for a generic pure state |<E>): 

t" = -($|r^|$) = 

TO 

However, for a free multiparticle system this is no longer true [3] because: 

i j i 

Also, for a single particle in a mixture of states the expectation value of r, i.e. tr(pr) is not vanishing unless the 
density operator is a scalar operator, i.e. commutes with all unitary representations of Lorcntz transformations. In 
fact, the operator exp[u; • J/T], that we have dealt with in this work, does not. 
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APPENDIX F - Polarization four-vector 



The results of Sect. V allow to determine the polarization four- vector of particles with momentum p and general 
spin S, in an ideal Boltzmann gas. This determination is alternative with respect to that presented in ref. [8]. The 
polarization four-vector is defined as: 



HAI W = ^MW^a) = ^((Wn) P (129) 

a 

where W is the Pauli-Lubanski four-vector: 

W = --€ J U PP T 

*' [J, 2 C HVpT° 1 

and p is the single-particle spin density operator, which, for our present case reads: 

P= ^yexp[-(e(p) + M g + o;-j)/T]P y 

being 

«(P) = V(p,<7|exp[-(£(p)+ M g + a;-j)/T]P y |p,a) = f d 3 x/(x,p) 
Thus, in view of (129), the polarization four-vector can be written in the inertial frame as: 

n = -^((I»p-- 

n(p) m 

n = -^«k» p x p + -^((j» P - (130) 

n(p) m n(p) m 

By using the results of eqs. (64), (65) and (77), we can readily find the various expression of the expectation values in 
the eq. (130), that is: 

«J» P = / d 3 x(xxp)/(x,p)+ / d 3 x/(x,p)4!^ 
Jv Jv X\t) 

«k» p = / d 3 x(Cp-ex).f(x,p) (131) 
Jv 



By plugging (131) into (130) we get: 



U = ^l-u: (132) 



which is the same result found in ref. [8]. 



